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Abstract 

Measurements of B~£ production and mass are performed with the decay mode 
B+ — > J/ipir + using 0.37 fb _1 of data collected in pp collisions at -^/s = 7 TeV by 
the LHCb experiment. The ratio of the production cross-section times branching 
fraction between the B+ — > J/ipTr + and the B + — > J/ipK + decays is measured to 
be (0.68 ± 0.10 (stat.) ± 0.03 (syst.) ± 0.05 (lifetime))% for B+ and B + mesons with 
transverse momenta j»t > 4 GeV/c and pseudorapidities 2.5 < r] < 4.5. The B+ 
mass is directly measured to be 6273.7 ± 1.3 (stat.) ± 1.6 (syst.) MeV/c 2 , and the 
measured mass difference with respect to the B + meson is M(B^) — M(B + ) = 
994.6 ± 1.3 (stat.) ± 0.6 (syst.) MeV/c 2 . 
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The B+ meson is unique in the Standard Model as it is the ground state of a family of mesons 
containing two different heavy flavour quarks. The inclusion of charge conjugate modes is 
implied throughout this Letter. At the 7 TeV LHC centre-of-mass energy, the most probable 
way to produce mesons is through the gg-fusion process, gg — > + b + c pQ. 

The production cross-section of the meson has been calculated by a complete order-a^ 
approach and using the fragmentation approach [Tj. It is predicted to be about 0.4 ]±b [2. 3J 
at = 7 TeV including contributions from excited states. This is one order of magnitude 
higher compared to that predicted at the Tevatron energy ^fs = 1.96 TeV. However, 
the theoretical predictions suffer from large uncertainties, and an accurate measurement 
of the B+ production cross-section is needed to guide experimental studies at the LHC. 
As is the case for heavy quarkonia, the mass of the B+ meson can be calculated by 
means of potential models and lattice QCD, where early predictions lay in the range from 
6.2-6.4 GeV/c 2 [I]. 

The B+ meson was first observed in the semileptonic decay mode B^ — > 
J/ip fi~)£ + X (£ = e,/i) by CDF [4]. The production cross-section times branch- 
ing fraction for this decay relative to that for B + — > J/ip K + was measured to be 
0.132 loioty (stat.) ± 0.031 (syst.) ±g;g|g (lifetime) for B+ and B + mesons with transverse 
momenta > 6 GeV/c and rapidities \y\ < 1. Measurements of the B+ mass by CDF [5] 
and DO [6] using the fully reconstructed decay B+ — >• J/0 (/-^aO 71 " 4 ^ ga ye M(B+) = 
6275.6 ± 2.9 (stat.) ± 2.5 (syst.) MeV/c 2 and M(B+) = 6300 ± 14 (stat.) ± 5 (syst.) MeV/c 2 , 
respectively. A more precise measurement of the B+ mass would allow for more stringent 
tests of predictions from potential models and lattice QCD calculations. 

In this Letter, we present a measurement of the ratio of the production cross-section 
times branching fraction of B^ — > J/ipir + relative to that for B + — > J/ipK + for B^ and 
B + mesons with transverse momenta pt > 4 GeV/c and pseudorapidities 2.5 < 77 < 4.5, 
and a measurement of the B+ mass, using 0.37 fb" 1 of data collected in pp collisions at 
yfs = 7 TeV by the LHCb experiment. The LHCb detector [7] is a single-arm forward 
spectrometer covering the pseudorapidity range 2 < r\ < 5, designed for the study of 
particles containing b or c quarks. The detector includes a high precision tracking system 
consisting of a silicon-strip vertex detector surrounding the pp interaction region, a large- 
area silicon-strip detector located upstream of a dipole magnet with a bending power 
of about 4 Tm, and three stations of silicon-strip detectors and straw drift-tubes placed 
downstream. The combined tracking system has a momentum resolution Ap/p that varies 
from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, and an impact parameter (IP) resolution 
of 20 |am for tracks with high transverse momentum. Charged hadrons are identified 
using two ring-imaging Cherenkov detectors. Photon, electron and hadron candidates are 
identified by a calorimeter system consisting of scintillating-pad and pre-shower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a 
muon system composed of alternating layers of iron and multiwire proportional chambers. 

The .B+ — > J/ip7T + and B + — > J/ipK + decay modes are topologically identical and are 
selected with requirements as similar as possible. Events are selected by a trigger system 
consisting of a hardware stage, based on information from the calorimeter and muon 
systems, followed by a software stage which applies a full event reconstruction. At the 
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hardware trigger stage, events are selected by requiring a single muon candidate or a pair 
of muon candidates with high transverse momenta. At the software trigger stage [HI E] , 
events are selected by requiring a pair of muon candidates with invariant mass within 
120 MeV/c 2 of the J/ip mass |10j . or a two- or three-track secondary vertex with a large 
track sum, a significant displacement from the primary interaction, and at least one 
track identified as a muon. 

At the offline selection stage, J/ip candidates are formed from pairs of oppositely 
charged tracks with transverse momenta pt > 0.9 GeV/c and identified as muons. The 
two muons are required to originate from a common vertex. Candidates with a dimuon 
invariant mass between 3.04GeV/c 2 and 3.14GeV/c 2 are combined with charged hadrons 
with px > 1-5 GeV/c to form the and B + meson candidates. The J/ip mass window is 
about seven times larger than the mass resolution. No particle identification is used in the 
selection of the hadrons. To improve the B+ and B + mass resolutions, the mass of the 
pair is constrained to the J/ip mass [TO] . The 6-hadron candidates are required to 
have pt > 4 GeV/c, decay time t > 0.25 ps and pseudorapidity in the range 2.5 < rj < 4.5. 
The fiducial region is chosen to be well inside the detector acceptance to have a reasonably 
flat efficiency over the phase space. To further suppress background to the B+ decay, the 
IP x 2 values of the J/ip and ir + candidates with respect to any primary vertex (PV) in 
the event are required to be larger than 4 and 25, respectively. The IP \ 2 is defined as 
the difference between the x 2 °f the PV reconstructed with and without the considered 
particle. The IP \ 2 °f the B+ candidates with respect to at least one PV in the event is 
required to be less than 25. 

The ratio of the production cross-section times branching fraction measured in this 
analysis is 

a(B+)B{B+^J/1>* + ) 



R 



c/u a{B+)B{B+ J/ipK- 

N(B+ -»■ J/^7T + ) 



(1) 



e t c ot N(B+^J/1>K+)' 

where a(B^) and a(B + ) are the inclusive production cross-sections of the B+ and 
B + mesons in pp collisions at y/s = 7 TeV, B(B^t —> J/ipn + ) and B{B + — > J/ipK + ) 
are the branching fractions of the reconstructed decay chains, N (B+ — > J/^7r + ) and 
N (B + ->■ J/ipK + ) are the yields of the ->■ J/ipTi + and B + -> J/ipK + signal decays, 
and e£ ot , e" ot are the total efficiencies, including geometrical acceptance, reconstruction, 
selection and trigger effects. 

The signal event yields are obtained from extended unbinned maximum likelihood fits to 
the invariant mass distributions of the reconstructed B+ and B + candidates in the interval 
6.15 < M(J/ipn + ) < 6.55GeV/c 2 for B+ candidates and 5.15 < M(J/ipK+) < 5.55GeV/c 2 
for B + candidates. The B+ — >■ J/ipir + signal mass shape is described by a double-sided 
Crystal Ball function [TT]. The power law behaviour toward low mass is due primarily 
to final state radiation (FSR) from the bachelor hadron, whereas the high mass tail is 
mainly due to FSR from the muons in combination with the J/ip mass constraint. The 
B + — > J/if) K + signal mass shape is described by the sum of two double-sided Crystal Ball 
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Figure 1: Invariant mass distributions of selected (a) — > J/vp7r + candidates and (b) B + — > 
J/ip K + candidates, used in the production measurement. The fits to the data are superimposed. 

functions that share the same mean but have different resolutions. From simulated decays, 
it is found that the tail parameters of the double-sided Crystal Ball function depend mildly 
on the mass resolution. This functional dependence is determined from simulation and 
included in the mass fit. The combinatorial background is described by an exponential 
function. Background to B + — > J/ipK + from the Cabibbo-suppressed decay B + — > J/i/jn + 
is included to improve the fit quality. The distribution is determined from the simulated 
events. The ratio of the number of B + — > J/ipir + decays to that of the signal is fixed 
to B(B+ J/^n + )/B(B + J/ipK + ) = 3.83% [12]. The Cabibbo-suppressed decay 
B+ — > J/ijj K + is neglected as a source of background to the B+ — > J/ip n + decay. The 
invariant mass distributions of the selected B^t — > J/ip7r + and B + — > J/ipK + candidates 
and the fits to the data are shown in Fig. [T] The numbers of signal events are 162 ± 18 for 
B+ ->• J/^7T+ and 56 243 ± 256 for B + ->- J/ipK + , as obtained from the fits. 

The efficiencies, including geometrical acceptance, reconstruction, selection and trigger 
effects are determined using simulated signal events. The production of the B + meson is 
simulated using Pythia 6.4 [13] with the configuration described in Ref. [33]. A dedicated 
generator BcVegPy [15] is used to simulate the B+ meson production. Decays of B+, 
B + and J/tp mesons are described by EvtGen [16] in which final state radiation is 
generated using Photos [17] . The decay products are traced through the detector by 
the Geant4 package [IB] as described in Ref. [TH] . As the efficiencies depend on and 
r), the efficiencies from the simulation are binned in these variables to avoid a bias. The 
signal yield in each bin is obtained from data by subtracting the background contribution 
using the sPlot technique [20] , where the signal and background mass shapes are assumed 
to be uncorrelated to pt and 77. The efficiency-corrected numbers of B+ — > J/ip tt + and 
B + — > J/ipK + signal decays are 2470 ±350 and 364 188 ±2270, respectively, corresponding 
to a ratio of R c / U = (0.68 ± 0.10)%, where the uncertainties are statistical only. 

The systematic uncertainties related to the determination of the signal yields and 
efficiencies are described in the following. Concerning the former, studies of simulated 
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events show that effects due to the fit model on the measured ratio R c / U can be as much as 
1%, which is taken as systematic uncertainty. The uncertainties from the contamination 
due to the Cabibbo-suppressed decays are found to be negligible. 

The uncertainties on the determination of the efficiencies are dominated by the knowl- 
edge of the lifetime, which has been measured by CDF |2T] and DO |22j to give 
t(B+) = 0.453 ± 0.041 ps [10]. The distributions of the B+ ->• J/ipn + simulated events 
have been reweighted after changing the B+ lifetime by one standard deviation around its 
mean value and the efficiencies are recomputed. The relative difference of 7.3% between 
the recomputed efficiencies and the nominal values is taken as a systematic uncertainty. 
Since the B + lifetime is known more precisely, the contribution to the uncertainty is 
neglected. 

The effects of the trigger requirements have been evaluated by only using the events 
triggered by the lifetime unbiased (di)muon lines, which is about 85% of the total number 
of events. Repeating the complete analysis, a ratio of R c / U = (0.65 ± 0.10)% is found, 
resulting in a systematic uncertainty of 4%. 

The tracking uncertainty includes two components. The first is the difference in track 
reconstruction efficiency between data and simulation, estimated with a tag and probe 
method |23j of J/ip — > decays, which is found to be negligible. The second is due 

to the 2% uncertainty on the effect from hadronic interactions assumed in the detector 
simulation. 

The uncertainty due to the choice of the (pt,v) binning is found to be negligible. 
Combining all systematic uncertainties in quadrature, we obtain R c / U = (0.68±0.10 (stat.)± 
0.03 (syst.) ± 0.05 (lifetime))% for B+ and B + mesons with transverse momenta pr > 4 
GeV/c and pseudorapidities 2.5 < r\ < 4.5. 

For the mass measurement, different selection criteria are applied. All events are 
used regardless of the trigger line. The fiducial region requirement is also removed. Only 
candidates with a good measured mass uncertainty (< 20 MeV/c 2 ) are used, and a loose 
particle identification requirement on the pion of the B^ — > J/ip7r + decay is introduced to 
remove the small contamination from B^ — > J /ipK + decays. 

The alignment of the tracking system and the calibration of the momentum scale are 
performed using a sample of J/ip —> decays in periods corresponding to different 

running conditions, as described in Ref. [21]. The validity of the calibrated momentum 
scale has been checked using samples of K® — > 7r + 7r~ and T — > decays. In all cases, 

the effect of the final state radiation, which cause the fitted masses to be underestimated, 
is taken into account. The difference between the correction factors determined using the 
J/ijj and T resonances, 0.06%, is taken as the systematic uncertainty. 

The .B+ mass is determined with an extended unbinned maximum likelihood fit to the 
invariant mass distribution of the selected B^ — > J/ip7r + candidates. The mass difference 
M(B^) — M(B + ) is obtained by fitting the invariant mass distributions of the selected 
B+ — > J/ipir + and B + — > J/il)K + candidates simultaneously. The fit model is the same 
as in the production cross-section ratio measurement. Figure [2] shows the invariant mass 
distribution for B+ — > J/ijj7r + . The B+ mass is determined to be 6273.0 ± 1.3 MeV/c 2 , 
with a resolution of 13.4 ± 1.1 MeV/c 2 , and the mass difference M(B^) — M(B + ) is 
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Figure 2: Invariant mass distribution of — > J/^7r + decays, used in the mass measure- 
ment. The fit to the data is superimposed. 

994.3 ± 1.3 MeV/c 2 . The uncertainties are statistical only. 

The mass measurement is affected by the systematic uncertainties due to the invariant 
mass model, momentum scale calibration, detector description and alignment. To evaluate 
the systematic uncertainty, the complete analysis, including the track fit and the momentum 
scale calibration when needed, is repeated. The parameters to which the mass measurement 
is sensitive are varied within their uncertainties. The changes in the central values of 
the masses obtained from the fits relative to the nominal results are then assigned as 
systematic uncertainties. 

Table [j] summarizes the systematic uncertainties assigned to the measured mass 
and mass difference AM = M(B+) — M(B + ). The main source is the uncertainty in the 
momentum scale calibration. After the calibration procedure a residual ±0.06% variation 
of the momentum scale remains as a function of the particle pseudorapidity rj. The impact 
of this variation is evaluated by parameterizing the momentum scale as a function of r\. The 
amount of material traversed by a particle in the tracking system is known to 10% accuracy, 
the magnitude of the energy loss correction in the reconstruction is therefore varied by 
10%. To quantify the effects due to the alignment uncertainty, the horizontal and vertical 
slopes of the tracks close to the interaction region, which are determined by measurements 
in the vertex detector, are changed by ±0.1%, corresponding to the estimated precision 
of the length scale along the beam axis (25]. To test the relative alignment of different 
sub-detectors, the analysis is repeated ignoring the hits of the tracking station between the 
vertex detector and the magnet. Other uncertainties arise from the signal and background 
line shapes. The bias due to the final state radiation is studied using a simulation based 
on Photos [17J. The mass returned by the fit model is found to be underestimated by 
0.7 ± 0.1 MeV/c 2 for the B+ meson, and by 0.4 ±0.1 MeV/c 2 for the B + meson. The 
mass and mass difference are corrected accordingly, and the uncertainties are propagated. 
The effects of the background shape are evaluated by using a constant or a first-order 
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Table 1: Systematic uncertainties (in MeV/c 2 ) of the mass and mass difference 
AM = M(B+) - M(B + ). 



Source of uncertainty 


M(B+) 


AM 


Mass fitting: 






- Signal model 


0.1 


0.1 


- Background model 


0.3 


0.2 


Momentum scale: 






- Average momentum scale 


1.4 


0.5 


- r] dependence 


0.3 


0.1 


Detector description: 






- Energy loss correction 


0.1 




Detector alignment: 






- Vertex detector (track slopes) 


0.1 




- Tracking stations 


0.6 


0.3 


Quadratic sum 


1.6 


0.6 



polynomial function instead of the nominal exponential function. The stability of the 
measured mass is studied by dividing the data samples according to the polarity of 
the spectrometer magnet and the pion charge. The measured B£ masses are consistent 
with the nominal result within the statistical uncertainties. 

In conclusion, using 0.37 fb _1 of data collected in pp collisions at y/s = 7 TeV by 
the LHCb experiment, the ratio of the production cross-section times branching fraction 
of B+ — > J/ipn + relative to that for B + — > J/ipK + is measured to be R c / U = (0.68 ± 
0.10 (stat.)±0.03 (syst.)±0.05 (lifetime))% for B+ and B + mesons with transverse momenta 
Pt > 4 GeV/c and pseudorapidities 2.5 < r\ < 4.5. Given the large theoretical uncertainties 
on both production and branching fractions of the B+ meson, more precise theoretical 
predictions are required to make a direct comparison with our result. The B£ mass is 
measured to be 6273.7 ± 1.3 (stat.) ± 1.6 (syst.) MeV/c 2 . The measured mass difference 
with respect to the B + meson is M(B+)-M(B + ) = 994.6 ±1.3 (stat.) ±0.6 (syst.) MeV/c 2 . 
Taking the world average B + mass |1Q| . we obtain the indirect measurement of the B+ 
mass M(£?+) = 6273.9 ± 1.3 (stat.) ± 0.6 (syst.) MeV/c 2 , which has a smaller systematic 
uncertainty. These results represent the most precise determinations of these quantities to 
date. 
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